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' Abstract 

(N . 

^jr^l When solving the equations of General Relativity in a symmetric sector, it is natural to 

^ ' consider the same symmetry for the geometry and stress-energy. This implies that for static 

I and isotropic spacetimes, the most general natural stress-energy tensor is a sum of a perfect 

fluid and a radial imperfect fluid component. In the special situations where the perfect fluid 
component vanishes or is a spacetime constant, the solutions to Einstein's equations can be 
thought of as modified Schwarzschild and Schwarzschild-de Sitter spaces. Exact solutions of 
this type are derived and it is shown that whereas deviations from the unmodified solutions 
can be made small, among the manifestations of the imperfect fluid component is a shift 
in angular momentum scaling for orbiting test-bodies at large radius. Based on this effect, 
the question of whether the imperfect fluid component can feasibly describe dark matter 
phenomenology is addressed. 
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^ ! 1 Introduction 

1/^ . The observed Universe is homogenous at very large scales but not at small ones where it contains 
numerous localized matter sources. Within the well-tested [1] theory of General Relativity [2], 
QO ' this means the geometry on large scales can be described by a metric of the Priedman-Robertson- 
Walker type while at smaller scales it must be approximated by metrics of the Schwarzschild 
I or Schwarzschild-de Sitter forms. These idealized metrics can be derived from the equations of 
motion of General Relativity by assuming different symmetries at the different scales. 

It may be argued that when solving the equations of motion of General Relativity in some 
symmetric sector, it is natural to consider the same symmetry for the geometry g^j/ and stress- 
ed I energy T^i, to which the geometry couples. Accordingly, when assuming spatial homogeneity 
for the geometry, it is natural to take T^i, to contain one free function associated with the time 
direction and one free function associated with the spatial directions. A stress-tensor of this 
kind is said to be in perfect fluid form and can be written as 

T^l = {p + p)U^U,+pgf,, (1) 

with a unit time- like vector, and p and p some functions called the energy density and 
pressure, respectively. In contrast, for static isotropic geometries, what is natural for T^ii is that 
it should have three independent functions - two associated with the time and radial directions, 
and one associated with the remaining two angular directions. A stress tensor of this type can 
be written as 

T^u=Tg + qV^K (2) 

where the first term on the right hand side is of perfect fluid form ([T]), and in the second term 
is a unit vector pointing in the radial direction and q is some function. Since the modification to 
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the perfect fluid form comes only in the radial components of T^,y, a stress-energy tensor of this 
type may be called radially imperfect. In other terminology it may be said to contain viscous 
shear. It is the most general stress-energy tensor compatible with staticity and isotropy [3]. 

In principle, the stress-energy tensor can be calculated in any theory of matter coupled to 
gravity given a description of the matter distribution. Recently, Cox, Mannheim, and O'Brien 
performed calculations of this kind assuming a matter component in the form of an incoherent 
collection of modes of a free (minimally-coupled) massless scalar field [4]. Reproducing earlier 
results [5], they showed that such a matter distribution puts the stress-energy tensor into the 
perfect fluid form when the background spacetime is Minkowski. Crucially, however, they also 
showed that the incoherent collection of fleld modes does not necessarily generate a perfect fluid 
stress-energy tensor if the background spacetime is nonhomogenous. In one analytic calculation 
involving a special static isotropic background, they put the stress-energy tensor into the form 
([2]) and explicitly wrote the functional form for p, p and q with q ^ 0. This corroborates the 
motivation for using the general form ([2]) for the stress-tensor source when it is coupled to a 
static and isotropic geometry. 

The calculations in were done in the spirit of field theory on curved spacetime. That 
is, the background spacetime was assumed to be fixed. The approach was useful because it 
allowed those authors to determine the field modes explicitly, sum over them, and thus obtain 
the stress-energy tensor. From the nonperturbative point of view, however, it is important to 
take into account backreaction, i.e. to obtain a consistent solution to Einstein equations with 
an imperfect fiuid source. The observations in [3] are a direct motivation for doing exactly that. 

The purpose of the present work is to study static and isotropic spacetimes coupled with 
the imperfect fluid stress-energy tensor ([2]). A similar task has been attempted in [3] but the 
issue is revisited here in more detail. Attention is focused on describing spacetimes in which 
the imperfect fluid component is nonvanishing and comparing them to known solutions, like 
Schwarzschild, in which this component is exactly zero. 

Section [2] deals with Einstein's equations for static and isotropic spacetimes coupled to certain 
fluids of the form ([2]). Although isotropy may in some systems also be thought to be an emergent 
property due to matter coupling [Bl [7] , in this work it is imposed strongly. Several solutions 
- modiflcations of the Schwarzschild, anti-de Sitter, and Schwarzschild-de Sitter spaces - are 
derived and their properties are discussed. Then, some phenomenological aspects of the modifled 
black-hole solution are studied in Sec. [3j Since one aspect of the phenomenology is similar to that 
attributed to dark-matter in galaxies. Sec. H] explores whether the imperfect fluid component 
can actually be considered a candidate description for dark matter. A brief summary of the 
results appears in SecHl 

2 Static Isotropic Solutions 

In units where the speed of light and the Newton constant are set to one, the equations of motion 
of General Relativity are 

G^u = S-nT^y (3) 

with G^y the Einstein tensor derived from a metric g^y, and T^y the stress-energy tensor of 
matter [2]. This section deals with solutions to these equations under the assumption that the 
geometry is static and isotropic (spherically-symmetric) and the stress tensor is in the radially- 
imperfect form ([2]). 
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With the assumed symmetries, the Hne element ds^ = g^y dx^ dx'^ associated with a metric 
Qfj^y in coordinates x^^ can be written as 

ds^ = -f dt^ + h-^ dr^ + dn'^ (4) 

so that the only unfixed components of the metric are gu = — and grr = h~^. The functions 
/ and h depend on the radial coordinate r only and = dO^ + sin^ 9 dcjp' is the line element 
of a unit sphere. With this choice of coordinates and parametrization, the nonzero components 
of the Einstein tensor occur only on the diagonal and are 

Gtt = ^{l-h- rh') (5a) 
Grr = ^{hf-f + 2rhf') (5b) 
Gee = ^{2h{f + rf") + h' {f + rf) ) (5c) 
G^^ = sm^eGee. (5d) 
With a radially imperfect fluid of the form ([2]), the equations of motion become 

(6a) 
(6b) 
(6c) 
(6d) 

where p, p and q are functions only of r, and factors of Svr have been absorbed on the right 
hand side. Since equations ()6d|) and ([5c|) are multiples of each other, there are actually only 
three independent equations. And since there are five unknown functions, the equations cannot 
be solved without specifying additional data. If the geometry, i.e. the functions / and h, is 
specified then the equations can be solved directly to yield the compatible p, p and q. In this 
way, any spacetime whatsoever can be obtained by some choices for the stress-energy tensor. 
But it is more constructive to specify reasonable sources, i.e. p, p and q, and subsequently solve 
for the compatible geometry. Below are a few cases characterized by idealized choices for the 
energy density p and pressure p. 

2.1 Exactly Vanishing Perfect Fluid 

The case of the vanishing perfect fluid is defined by 

P = P = (7) 

everywhere. In this case, equation (j6a|) can first be solved for h. Next, the solution for h can be 
inserted into ([6cj) to yield a linear equation for /. The results are 

(8) 
(9) 
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= (p + q) h~ 
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where a, b, c, d are constants, a and c are dimensionless, while b and d have units of length. (The 
constant c should not be confused with the speed of light which is set to one throughout by a 
choice of units.) These functions / and h determine the metric. When they are inserted into 
(|6bp. they determine the imperfect fluid component to be 



If the imperfect fluid component is required to vanish, then the constant c must be set to 
zero. When in (l8|) and / in ([9]) with c = are inserted into they yield the Schwarzschild 
line-element with an additional numeric constant in front of the dt"^ component. Because 
of reparametrization invariance, this constant can be set to unity without loss of generality. 
Thus the standard black hole solution, with a horizon at the Schwarzschild radius r = b and a 
singularity at r = 0, is recovered. 

If c 7^ 0, the coefficient of dr'^ in the line element remains identical as in the Schwarzschild 
spacetime. In this sense the quantity b may still be called the Schwarzschild radius. The 
coefficient of dt'^ in the line element, however, is for c 7^ very different than before. This can 
be seen in Fig. [1] where it is plotted for a few combinations of parameters. A number of new 
properties are worth pointing out. 

First, the functions have a root at a finite radius that is different from the usual r = b. 
Because of the dependence of q on f, the imperfect fluid component diverges at those roots. 
The curvature scalar R^uXaR^'^^'^ also diverges there and so these points (actually, surfaces of 
spheres) are likely unphysical features of the solutions. Their location moves close to r = 6 as 
c approaches zero through positive values, and away from r = b as c approaches zero through 
negative values. If the purpose of the spacetime is to describe behavior outside the black hole, 
positive c should probably be preferred. 

Second, the solutions are only valid for r > b and it is not clear whether they can be extended 
to r < 6. Moreover, the quantity is not zero at r = fe so even if the solution is continued 
to smaller values of r the spacetime signature may change. Although r = 6 is always "hidden" 
behind the singularity described above, this is nonetheless an unphysical feature that suggests 
that the spacetime should only be studied far from r = b. 

Third, for large r, the functions diverge rather than approach a constant and so they do not 
describe asymptotically flat spacetimes. This occurs even though the nonvanishing component of 
the stress-energy tensor, proportional to q, falls to zero as r — > 00. This is yet another undesirable 
feature because it implies that the spacetime, even though it is a modification of Schwarzschild 



4 



space, does not have a limit in which it reproduces Newtonian physics. 

As a special case, it is possible to consider the solution with 6 = 0. In the Schwarzschild 
solution, this implies removing the black hole from the spacetime leaving only flat space. In the 
case with c / 0, setting 6 = gives 

/4r" 

a-c + clogW— . (II) 

While this is valid for all r > 0, the other two features of the modified black hole spacetimes 
mentioned above persist. Namely, the function has a root at some finite value of r where the 
stress-energy tensor component q as well as the scalar curvature R^j^yXaR^^^" diverge, and it 
does not asymptote to a constant for large r. Thus, the resulting spacetime represents quite a 
drastic departure from Minkowski spacetime and should be interpreted similarly as above. 

2.2 Cosmological Constant 

The case of the cosmological constant is defined by a perfect fluid with 

P = -P = K (12) 

where A is a nonzero spacetime constant. The strategy for obtaining solutions in this case is 
the same as before. Equation (j6ap can first be solved for h to give 

6 Ar2 

h = l 13 

r 6 



with 6 a constant. Then, this can be substituted into (j6cj) to give a linear equation for /, 

2r2 (36 - 3r + Ar^) /" + r (36 - 6r + 4Ar3) /' - (36 + 4Ar3) / = 0. (14) 



In the special case 6 = and A < 0, the solution to ()14p can be found analytically and then 
inserted into (|6b|) to give q. The result is 

arctanh ( a/ 1 - ^ 1 | , (15) 





(16) 

with a and c constants. 

The effect of the imperfect fiuid can again be removed by setting c = 0. The resultant 
function / with c = and a = 1 set by convention, when inserted into the line-element ansatz, 
together with (jl3p with 6 = 0, gives anti-de Sitter spacetime. This is actually a homogenous, 
rather than only isotropic, space because 6 = 0. 

If c / 0, the solution describes a departure from anti-de Sitter. Figure [2] shows the quantity 
= —gtt for a few choices of parameters. For c > 0, the curves have a root at some finite radius 
where q and the curvature invariant R^vXaR^^'^" diverge. Curves with c < are nonsingular 
everywhere. For both signs, the curves significantly differ from the anti-de Sitter curve (c = 0) 
near r = but track it closely for r ^ oo. The modified spacetimes are thus asymptotically 
anti-de Sitter. 
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For A > 0, it can be checked that the function 



is a solution to and corresponds to g = 0. It describes Schwarzschild-de Sitter spacetime 
a b ^ and a de Sitter spacetime if 6 = 0. A feature of these spaces is that they contain 
a cosmological horizon at a radius rc for which fsdsii^c) = 0. The metric ansatz in static 
coordinates suffers from a coordinate singularity at this radius as both gu and grr change sign and 
grr diverges. Because of this solutions must be specified separately for regions inside and outside 
the cosmological horizon. While this is not a serious hurdle when g = 0, it does complicate 
matters when q ^ 0. For example, (jl5|) and p6|) can be used for A > and 6 = to describe 
solutions outside the horizon but they cannot be used in the region inside the horizon because 

is not real-valued everywhere there for any choice of parameters a and c. A solution inside 
the horizon does exist - it can be computed numerically given some boundary data - but it is 
not described by (fTSl) and ([TB]) . 

At this stage, it is worth noting that the function q in (jlOp and (jl6p has the same form in 
two inequivalent situations. It seems reasonable therefore to guess that this form is also valid 
in the more general case for A > and b ^ 0. With this guess, the general solution for / can 
be obtained by solving equation (I6b|) and then verified by substitution into (pcj) and (fT4|) . The 
guess works and the result is 



b Ar^ 

/ = ayi- -- — -c/e,A,& (18) 

Q = 4? (19) 

where fc,A,b is a function that can be written analytically but has a complicated form that does 
not provide insightful information. 

Fig. [3] shows = —gtt obtained using p8|) with various parameters. All curves are plotted 
in the domain of r between the two horizons of the Schwarzschild-de Sitter solution. Some of the 
features that distinguish them from the Schwarzschild-de Sitter case are qualitatively similar to 
those already found in the modified Schwarzschild spacetimes (A = 0). 

First, all curves with c 7^ have zeros at radii that are different from the black- hole and 
cosmological horizons in the unmodified spacetime. The curvature invariant R^vXaR^'^^^ diverge 
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at those zeros. Their location moves close to the black-hole horizon as c approaches zero through 
positive values, and move close to the cosmological horizon as c approaches zero through negative 
values. 

Second, for c 7^ 0, the curves do not fall to zero at the black hole or cosmological horizons. 
This poses problems with signature change. 

Third, the maxima of the curves with c 7^ are shifted relative to the c = case. For 
c > 0, the maximum occurs closer to the cosmological horizon. For c < 0, it occurs closer to the 
black-hole horizon. 



2.3 Almost-Constant Perfect Fluid 

This case summarizes and elaborates on the discussion in [3]. It involves a more qualitative 
description of the effect of the imperfect fluid component on geometry when the assumption 
that p and p be strict spacetime constants is relaxed. The profiles for p, p and q are assumed to 
satisfy 

^-{f* till 

_ r q{r) for r < R, 

~ \ for r > i?g ^^^^> 

with appropriate continuity conditions at r = Rp,Rp,Rq. The boundary radii define compact 
regions where the corresponding stress-energy tensor components differ from constants. A cos- 
mological constant is allowed for large r in the profiles for p and p but not for q. 

When (j20ap is inserted into (j6ap . the resulting equation can be solved in the region r > Rp 
to give the solution 

b Ar^ , , 

h=l for r>Rp. (21) 

r 3 

This can then be inserted into ()6bp and ([6c|) . The results that ensue depend on the relative sizes 
of Rp, Rp and Rg. There are six possible orderings and these can be treated exhaustively. 

In general, the strategy for obtaining exact solutions is to solve (j6c]) for r > Rp and r > Rp 
where the details of the p and p profiles are not important. The equation, and therefore the 
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Figure 4: Isotropic spacctimc witli p, p and q profiles with compact support. 



solution, in this region are equivalent to that studied in the previous sections. When the solution 
is inserted into ()6bp . the necessary form for q is obtained and it is not a function with compact 
support as required by (j20cp . To be consistent, then, it is necessary to require that q = in 
that region. Thus, the exterior solution becomes 

/ = aVh for r>Rp, r> Rp. (22) 

The pathologies associated with large radii described for the idealized solutions in Sec. 12.11 and 
12.21 are thus removed. 

An interesting case is when Rp = Rq > Rp, illustrated in Fig. HI As in the general case, 
the solution in the exterior region r > Rp reduces to ()22p . However, in this case there is a shell 
Rp < r < Rp where 

/ / aVh. (23) 

This general feature may also be argued to hold in regions r < Rp where (? 7^ [8]. 

The precise form of / in the shell Rp < r < Rp depends on the pressure profile p. Under 
certain conditions, however, it may be reasonable to approximate / using ([9]) or (jlSp when p 
does not differ significantly from a constant. This can only be done in a region from r = Rp to 
some intermediate radius shown by the dashed line in Fig. [H Beyond this intermediate radius, 
the pressure must differ sufficiently from a constant in order to ensure that / goes to (I22p for 
r > Rp. 

3 Phenomenology 

The most interesting among the solutions with imperfect fluids is the modified black hole space- 
time with compact source of Sec. 12.31 However, that solution cannot be written explicitly 
without assuming profiles for p, p and q. In contrast, the modified black hole spacetime of 
Sec. 12.11 can be written down analytically in terms of only two free parameters. Moreover, as 
argued in Sec. 12.31 it may be a reasonable approximation to the more physical spacetime in a 
region close to the matter distribution with energy density p and away from r = 00 (or r = 
if A > 0). So it may be instructive to consider that idealized spacetime as a concise model of 
the possible effects of a radial imperfect fluid component on geometry. The model does have 
limitations, however. In particular, it does not easily allow to estimate effects due to multiple 
source. Bearing such limitations in mind, this section aims to compare the phenomenology of 
the modified spacetime, following [2], with the Schwarzschild black hole. 
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The line element of the modified Schwarzschild solution, repeated here for convenience, is 

ds^ = -f dt^ + h-^ dr^ + dn^ (24) 

with 

/ = alWl--)-c(l-Wl--log(^i— -] 1 (25) 




o'^' + Veff = ^E' (28) 



h = l--. (26) 
r 

It reduces to that of Schwarzschild spacetime for a = 1 and c = and describes modification 
thereof when c ^ 0. In the analysis below, all analytic results are shown with a explicitly present 
but all numerical quantities are evaluated using a convention setting a = 1. 

3.1 Orbital motion 

A unit four- velocity vector = dx^/dr, where r is the proper time, should obey the normal- 
ization condition —1 = u^u^. Assuming motion in the plane 6 = '7r/2, this condition gives 

-1 = - ff + h~^r^ + r^^'^ (27) 

with the dot denoting d/dr. Because {d/dtY and [d/dcf))^ are Killing fields of the static isotropic 
spacetime, the quantities E = f'^i and L = r"^ (p, called energy and angular momentum, respec- 
tively, are constants of motion. Writing (j27p in terms of E and L yields 

1.0 1 
— 7 

2 

with 

Veff = ^ (l + 5) + {^-^)- (29) 
This is an equation for a test particle moving in one dimension in an effective radial potential 

Veff. 

A test particles moving in a circular orbit is at the minimum of the potential and has 
dVgff/dr = 0. Since it does not move radially, it also has r = so that Veff = E'^/2. Together, 
these two conditions can be used to solve for L and E in terms of the orbit radius. The general 
result 

E^ = — ^ ~ a2 + Gc log- (30) 
f -rf d 

is independent of the function h. The last approximations are obtained after expanding for large 
r and small c and ignoring all but the leading terms that are either proportional to c or are 
independent of it. 

The general result reproduces the Schwarzschild scalings when c = 0. For c 7^ 0, however, 
the correction terms can dominate for r sufficiently large. In particular, the scaling of angular 
momentum squared, L^, shifts from linear to quadratic near 



9 



Since angular momentum L is the product of radius r and velocity f , a scaling oc 
implies that the velocity of orbiting particles is independent of distance. Interestingly, this 
behavior is actually observed within galaxies where it is interpreted either as evidence for the 
existence of dark matter or as a motivation for modifying the dynamics of general relativity. 
Phenomenologically, the observed behavior is well-described by the framework of Modified New- 
tonian Dynamics (MOND) [9l [10] which postulates that orbiting bodies deviate from standard 
Newtonian behavior and start orbiting at a constant velocity when their acceleration falls below 
the threshold 10~^^m~^. (And by what seems a coincidence, the acceleration threshold is within 
an order of magnitude of the root of the observed cosmological constant A.) Setting this equal 
to the Newtonian expression implies that the threshold radius in the MOND description is 

r^MOND = (3-10^3 ml/2) (33) 

If the derived threshold (j32p is to account for galaxy rotation curves, then, assuming a fixed 
parameter a, the parameter c should depend on b as 

c = 7\/6 (34) 

with 7 some quantity with units of mT^I"^ . Substituting this relation into (j32p and setting that 
equal to ([33]) gives 

7A,f = 3-10-1^ m~i/2_ (35) 

The subscript denotes this value is obtained from the MOND phenomenological description. 

If the derived threshold ([32p is said not to be responsible for galaxy rotation curves, then 
a bound on c can be formulated from each galaxy's data. It is convenient to summarize those 
bounds using the parameter 7 from ([34p . Since the size of a galaxy is typically an order of 
magnitude larger than (|33p . the bound on 7 would in this case be about an order of magnitude 
lower than ([35|) . 

Within the solar system, objects orbit the Sun according to the unmodified scaling at least 
up to radii 50 times the Earth's orbit radius. Using h as the Schwarzschild radius of the Sun, 
these observations yield the bounds 

\cs\ < 10^1°, I75I < lO^i^m^i/^ (3g) 

The subscripts denote that the bound refers to an observation wherein the central massive object 
is the Sun and the bound on 7 is presented assuming the scaling ([M[) . 



3.2 Perihelion Precession 

If a test particle is displaced from a circular orbit, it oscillates around the equilibrium radius 
with frequency ujr given by = d?Veff/dr'^ with the right-hand-side evaluated at the orbit 
radius. Its angular frequency is = L/r^ evaluated at the orbit radius. When these are not 
equal, the object precesses with frequency = tj^—ujr- Using the modified black-hole spacetime 
and evaluating cOp for large r and small c yields 

~ ^ c ■= (37) 

The precession frequency for an elliptic but close to circular orbit differs from this result by a 
factor of order one [2]. 
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The precession of the planet Mercury around the Sun is consistent with the c = form of 
ujp with an accuracy of 10~^ pT]. This leads to 

\cs\ < 10"^^ I75I < 10-2^m-^/2_ (3g) 

Again, the subscripts denote that the central source is the sun. The bound on 75 is obtained 
using pill . 



3.3 Redshift 

In the static coordinates of (j24p . = (1,0,0,0)"^ is a time-like Killing vector. Its magnitude is 
{—(.fi^,^)^^"^ = f and so a unit four-velocity vector parallel to can be written as = 
The frequency a; of a light ray with momentum vector recorded by an observer with four- 
velocity is u; = k^u^ = fc/^C^- If two observers with unit four-velocity located at two 
different locations labeled by radii ri and r2 measure frequencies iOi and 0^2 for the same light 
ray, the ratio of their results is then 

^ = f4. (39) 

Without loss of generality, the coordinate r2 can be replaced hy r2 = ri{l + 6) for some 5 > —1. 
It is convenient to also define a quantity 

Y = '^^l^ = l-^ (40) 
L02 0J2 

that measures the relative change in the measured frequencies. This can be computed in the 
limit of large ri , small c and small 5 to yield 

Y ~ (41) 

after omitting terms of order c^, r~^, 5^, c/r, and higher. 

In Earth-based redshift experiments, the accuracy with which observations match (|4ip with 
c = is 7-10~5 p:2]. This gives 

IcijI < 10"l^ |7s| < 10~^2^~l/^ (42) 

where the subscripts indicate the result refers to an Earth-based observation. 



4 Discussion 



In a generic theory of the universe, the matter field content can be split into standard-model 
fields ■0, a scalar field 99, and other fields x thought to describe dark matter. Each of these 
fields can in principle give rise to a stress-energy tensor with an imperfect fluid component 
q. Sec. [3] shows that the effects of this component in central-mass situations can be made 
as small as necessary by changing the parameter c. The argument from Sec. 12.31 also implies 
that the imperfect fluid effects can be totally removed at large scales while producing effects 
near conventional matter sources that are similar to those usually associated with dark-matter. 
This section thus aims to combine theoretical ideas and phenomenological constraints in order 
to understand to what extent and under what conditions the imperfect fluid component can 
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account for dark-matter effects and thereby dispense with the need to postulate the existence 
of separate dark-matter fields x- 

For simphcity, attention can be concentrated on the scalar field (p. A large class of theories 
involving such a field have an action of the form 




where V{'^) is a potential that does not depend on derivatives of ip. The stress-energy tensor, 
obtained by variation with respect to the metric, is 

Tf,y = d^ifdyif - ^gf,ug°''^da(pdpf - g^,vV{ip). (44) 

For a field configuration with d^^p = 0, the tensor becomes T^y = —g^yV{ip) and thus describes 
a cosmological constant with A = V{ip). On small scales, it is plausible for the field ip to respond 
to the curvature of spacetime around massive objects, break the condition d^(p = in certain 
regions, and thus lead the stress-energy tensor to acquire more structure. In static and isotropic 
regions, the tensor will take the general form ([2]) with the energy density p, pressure p, and 
imperfect fluid component q related to the derivatives of the field. In this way, the component 
q acting like a dark field can become a companion of a conventional energy density p [14^ I15j . 

As shown in Sec. [21 the imperfect fluid component g in a static isotropic spacetime takes the 
form q ~ c/r'^f. The quantity c that determines its magnitude is a real dimensionless number 
and it is not unreasonable for it to actually represent a ratio of physically relevant scales. The 
dimensionful quantities in the general central source setup of Sec. 12.31 are A, p, p, Rp, Rp, and 
Rq (henceforth the quantities p and p are taken to represent averages of energy density and 
pressure over the relevant regions of spacetime) . The quantity b, the Schwarzschild radius of the 
energy density distribution, is a derived quantity proportional to pRp. In the situation depicted 
in Fig. m the radii obey Rp ~ Rq > Rq and it is reasonable to assume that Rp, Rq and p are 
somehow related to A. The independent dimensionful quantities are thus reduced to A, and two 
of p, Rp, and b. 

Other dimensionful parameters may enter from the specific theory describing the field ip. 
In a quantum theory, for example, a new dimensionful quantity is the Planck mass Mp. In 
an interacting theory, dimensionful parameters may also appear from coupling constants in 
potentials involving ip. For simplicity, however, the following arguments assume that these 
particle-physics parameters are not important. 

The most general dimensionless quantity that can be composed of the relevant parameters 
can then be written as 

= 6"" A"'^ p""" = p A''^ (45) 

where a and (5 are real numbers either of order one or with a physical interpretation, and the 
exponents obey 

= Uh — 2nA — 2np = Ub — 2?ia + n/j. (46) 

The discussion of orbital motion in Sec. [3] suggests that the scaling should be c = jVb, which 
corresponds to setting rih = 1/2. 

The value of 7 appears to be (|35p in galaxies but the bound (j38p from Mercury perihelion 
precession indicates that it must be significantly smaller within the solar system. Thus, if 7 is 
a universal constant, the perihelion precession observation rules out using the imperfect fiuid 
effect on orbital motion as a description for dark matter leading to flat galaxy rotation curves. 
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If 7 is not a universal constant, then the powers Up or n/j in (j45p should be different from 
zero. If TiA > 1/4 and < 0, the effective value of c can be much different in galaxies and 
within the solar system. For example, for nA = 1, c becomes 

c = ab'/^Ap-^/'' = (3b'/^AR'/\ (47) 

If p ~ A and a ~ 1, then c ~ and the boundary between linear and quadratic angular 

momentum scaling occurs at a radius close to (I33p . A similar boundary radius can be obtained 
by taking a larger and more realistic value for p for a galaxy and multiplying by a larger number 
a. For a dense system, p ^ A, the value of c becomes much smaller and the transition radius 
()32l) correspondingly larger. In this way, the imperfect fluid model with such scaling can be 
made to account for galaxy rotation curves and also be consistent with all the observational 
bounds described in Sec. El 

In principle, the correct scaling should be computed from some theory using (|44p . Such a 
computation would be desirable also because it could shed light on how multiple sources affect 
the strength of the imperfect fluid in realistic mass distributions. However, this requires solving 
for and manipulating the modes of the field (p in the curved background described in Fig. HJ 
and is not straight-forward. Without the specific theory, the rather special form of the scaling 
must be regarded only phenomenological fit. 

If the interpretation of dark matter effects in terms of an imperfect fluid related to (p were 
to hold, it would have a number of consequences. First, it would imply that direct experimental 
searches for dark matter particles should give results consistent with properties of the field 
ip responsible for the background cosmological constant. Second, the required scaling for the 
parameter c would imply that the threshold radius should depend on the density or size of a 
galaxy. This effect may be testable once the imperfect fluid is better understood. Third, it would 
imply that the pressure component of the field responsible for q cannot be constant everywhere. 
If it were, the imperfect fluid phenomenology would not be localized around the central source 
and would be severely inconsistent with large-scale properties of the observed universe. 

Other phenomena such as gravitational lensing and structure formation in cosmology [101 
[T3| \TB[ [TT] may provide more constraints when analyzed in the context of the imperfect fluid. 
Distinguishing between a dark-matter field source and an imperfect fluid component should be 
possible because because the two have slightly different signatures: whereas the former affects 
both the grr and gu components of the metric, the imperfect fluid affects primarily gu- Related 
issues have been studied in the context of alternative theories of gravity [lOl [13] . 

5 Conclusion 

Static and isotropic geometries can be consistently coupled to a stress-energy tensor composed 
of a perfect fluid component plus a radial imperfect fluid component. In cases where the perfect 
fluid component vanishes exactly or describes a cosmological constant, the solutions to the 
Einstein's equations are modifications of black hole spacetimes. One of such analytic solutions, 
described in Sec. 12.11 can be thought of as a modification to the Schwarzschild black-hole. Its 
line element differs from that of Schwarzschild space in the dt'^ coefficient, which in addition to 
the Schwarzschild radius b also contains two new parameters: one of these, c, is dimensionless; 
the other, d, has dimensions of length or mass. The spacetime contains a curvature singularity 
at r > 6 and is not asymptotically flat. Both these properties are pathologies that can be traced 
to an assumption that the pressure of the perfect fluid component of the stress-energy tensor 
vanishes exactly, and as discussed in Sec. 12.31 may be removed by relaxing this assumption. 
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Despite its limitations, the modified Schwarzschild spacetime can nonetheless be useful in 
describing how a radial imperfect fluid component can affect geometry. If the modified black- 
hole spacetime is treated as a model for the Earth's or the Sun's gravitational field, then, as 
Sec. [3] shows, various solar system observations constrain the parameter c. 

Among the deviations from Schwarzschild phenomenology is a shift in scaling of angular 
momentum L of an orbiting test body from ~ r to ~ at large r. Because this kind of 
behavior is an observed feature in many galaxies and one of the motivations for studies involving 
dark matter or modified gravitational dynamics. Sees. [3] and H] discuss conditions under which 
the imperfect fluid model can be used as a description for dark matter while still being consistent 
with solar system constraints. The conclusion is that for this to be possible the parameter c 
must scale with various other quantities in a rather special way. The scaling makes the setup 
somewhat contrived but also falsifiable. 

That a new term in the stress-energy tensor can produce new phenomenology should not in 
itself be surprising. After all, it is an additional source in the equations of motion and thus affects 
geometry. The description of dark-matter effects in terms of an imperfect fluid component, 
however, may be attractive in some regards. For example, if an imperfect fluid component of 
the correct magnitude were indeed produced by a scalar field 93 associated with the cosmological 
constant, it may help understand dark-matter effects without the need to introduce otherwise 
undetected fields or to modify gravitational interactions. From this perspective, the results from 
Sec. |4] suggest that the imperfect fluid component may at least be one of several complementary 
aspects to the wider subject of dark matter. Alternatively, that discussion can be viewed as a 
set of conditions that an effective theory of a dark-energy field ip should satisfy for it not to be 
at odds with solar-system and galaxy phenomenology. 
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